Highly interactive kink solutions 
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In this work we present a new class of real scalar field models admitting strongly interactive kink 
solutions. Instead of the usual exponential asymptotic behavior these topological solutions exhibit 
a power-law one. We investigate the interaction force between a pair of kink/anti-kink solutions 
both analytically and numerically, by integrating the time dependent field equations of the model. 
Furthermore, working within the first-order framework, we analyze the linear stability of these 
solutions. The stability analysis leads to Schodinger-like equations with potentials which, despite 
admitting no bound states, lead to strong resonance peaks. We argue that these properties are 
important for some possible physical applications. 

PACS numbers: 11.27.+d, ll.10.Lm 



I. INTRODUCTION 

Topological structures associated with solitary wave 
solutions of non-linear field theories are of great interest 
in several areas of physics PHI]- Particular topological 
(3, l)-dimcnsional kink-like defects, associated with in- 
terfaces between distinct order parameter regions, have 
been studied in detail both in non-relativistic and rela- 
tivistic contexts. Phase-field models have been success- 
fully applied for describing non-relativistic interface dy- 
namics across many fields in condensed matter and ma- 
terial science systems 043 ■ The dynamics of relativistic 
interfaces has also extensively been investigated, mainly 
in the context of cosmology, where domain wall networks 
that should have been formed in primordial phase transi- 
tions are considered, along with their cosmological impli- 
cations [l^-O- Recently a unified theoretical paradigm 
for interface dynamics, which includes both relativistic 
and non-relativistic systems in a unified framework, has 
been proposed in 14|. Moreover, in (4, l)-dimensional 



at the minima. As it is known, the second-order deriva- 
tive of the potential is related with the mass scale of the 
quantum meson, associated with the frequency of oscil- 
lation of the field at the minimum of the potential. A 
linear perturbation 77 of a vacuum solution, <p = v + rj(x), 
obeys the fluctuation equation Or] + m 2 rj = 0. The per- 
turbation solution, 77 = cos(u;f — kx), oscillates with an 
energy u 2 = k 2 + m 2 where the meson mass is given by 

ml = V^ 2 \v) (the notation V^ n) = d n V '{<f>) / 'd<f> n , for 
integer n is being used). In models admitting kink solu- 



tions 
is given by 
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^ 0, the asymptotic solution 



(x) ~ v + (x) where 



dx 2 



(1) 



braneworld theories, one can construct thick Minkowski 
branes which splits in a warped bulk spacetime due to a 
fist-order transition in a warped bulk [151 ] . 

In the present work we focus our attention on kink- 
type solutions, as topological static solutions of models 
described by a single real scalar field in bidimensional 
spacetime. Relativistic scalar fields models, admitting a 
large variety of kink- type solutions, have been studied 
in a diversity of scenarios [l6l - |23j . Models allowing dif- 
ferent kink/anti-kink interaction behavior [24l. |25l l27l - [29| 
have been motivated both theoretically and experimen- 
tally. Peculiar kink dynamics have also been studied in 
the context of models with modified kinetic term [30l - [35j 

We introduce a new model with the particular feature 
of having a null second-order derivative of the potential 



and 4>* ( x ) has an exponential behavior 0* (x) oc e m " x 
when x — > 00. When V} n '(v) = 0, all the linear per- 
turbation modes are massless and the solutions asymp- 
totic behavior is related to a higher order of the potential 
derivative at the minima. It is possible to show that for 
a model, in which the first n order derivatives of the po- 
tential vanish at the minimum, the equation of motion 
for the asymptotic kink solution is given by 



dx 2 



a n +i4>* 



(2) 



where a n = vl (v) , allowing a power-law behavior of 
model's asymptotic solutions, 



An 



where the constant A n depends on n and a r , 



(3) 
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In this paper we introduce the smooth minima poten- 
tial model corresponding to n = 3. The power-law behav- 
ior of this model asymptotic solution, and correspond- 
ing energy density, originates a non-standard kink/anti- 
kink interaction pattern, that is analyzed in detail in this 
work. We also investigate the generalized n-th depen- 
dent class of models finding for the fluctuations narrow 
resonance peaks whose lifetime gets reduced with the in- 
creasing of n. 

The paper is organized in the following way: In section 
[II] we explore the framework of the scalar fields theories 
that support kink solutions. In section 11111 we introduce 
the model and discuss its specific properties. In section 
IIVI we investigate the linear stability of the model so- 
lutions. In section [V] we present the generalized class 
for these models and we discuss resonance effects. We 
end the work in section IVII presenting some comments, 
conclusions and perspectives. 



II. THE FRAMEWORK 



Let us consider the standard Lagrangian density 



£ = i^W-V(0), 



(4) 



where <fi is a real scalar field and V(<j)) is the poten- 
tial which specifies the model under consideration, x M = 
(t,x) and x^ = (t, — x). We consider that the potential 
engenders a set of critical points, = {<j>i,(j>2, ■ ■ ■ , 4> n }, 
such that (m) = and Vfa) = for i = 1, 2, . . . , n. 
The equation of motion following from the above model 

is 



dt 2 dx 2 



vr = o. 



(5) 



The corresponding components of the energy- 
momentum tensor are given by 



(6b) 



For static solutions = 4>(x), we can write the equa- 
tion of motion as 



— = vi x) 

dx 2 



Integrating the equation of motion we obtain 



2 £ = vw+c 



(7) 



(8) 



where C is an integration constant, identified as the solu- 
tion pressure. The topological solutions necessarily obey 
the boundary conditions 



d<b 

lim -j— — >■ 

x— y±oo dx 



(9) 



in order to ensure finiteness of the gradient portion of 
the energy. These solutions are stressless, guaranteeing 
C = and V"(u,) = for the vi set of potential minima. 
For C = and V(4>) > 0, Eq. © is then given by 



The first (second) equation supports the monotonically 
increasing (decreasing) solution connecting two adjacent 
minima, called kink (anti-kink). For simplicity we will 
restrict ourselves to the explicit kink-type solution. The 
anti-kink solution may be obtained by the transformation 

4>{x) = 4>{-x). 

The superpotential function W(4>) may be defined as 



W(4>) = J y/2V(4>)d<l>, (11) 
allowing to express the energy of the solution as 

E : / p(x) = \W(4>(p6)) - W(^(-oo))| (12) 



without knowing the explicit form of the solution. 

The first order equations (TTU)l . for the kink solution, 
may then be written as [H, l37j 



dx 



= W. 



(i) 



(13) 



A model extensively studied in various scenarios 
[H, HH is the well-known double- well potential given by 



vM) = \{b 2 -tf)\ 



(14) 



where b and A are positive parameters. (±6) = 4A6 2 
The superpotential function for this model is then given 
by 



leading to the static solution 

4> s (x) — 6tanh(o \/Xx), 
and with an energy density given by 

p s {x) = A6 4 sech 4 (6VAa;). 
The corresponding total energy of the solution is 

ti x = . 



(15) 



(16) 



(17) 



(18) 
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The thickness of solution (fTB)) is defined by the following 
condition 5 S = 2(b\f\)~ 1 , such that 



4>s(oo) 



tanh(l) w 0.762. 



(19) 



The solution is fully characterized by its energy and 
thickness. 



III. THE MODEL 



We introduce the following model 



v(4>) 



A 



/ 2 13 



(20) 



where a and A are positive parameters. The model ex- 
hibits a Z2 symmetry tf> — > — cf>, with a local maximum 
point in <fi — 0, for V(0) = Aa 6 /2, and two minima points 
at (f> = ±a. It is possible to observe in Fig. [TJ that 
the two potential minima, with the particular character- 
istic of having a null potential second-order derivative 
V} (±a) = 0, are smoother than in (IT4l) model. 




FIG. 1: Potential profile of model (|14|l (solid line) and model 
pD| (dashed line) for a = b = X = A = 1. 



The model superpotential function is given by 



W(4>): 



with 



VXa 4 



/(0) + 3 arcsin — 
—f{4>) + 3arccosh 



101 



3tt 

T 



/(0) = ^>(5a 2 -20 2 )vV_^|. 



< —a 

|0| < a 

(f) > a 
(21) 

(22) 



The values of the superpotential function at the minima 
are given by W{±a) = 3y/\a 4 /16. 



A similar model was presented in Ref [40j, as a defor- 
mation of (j) 4 model (JH]), however in that case the po- 
tential function is positive defined and the correspond- 
ing kink and anti-kink solutions connect two inflection 
points. 

The equation of motion of model (I2U1) is given by 



3\0(a 2 - (f> 2 )\a 2 - (f> 2 \ =0, 



with the corresponding static solutions, 
ing the following equation 



^ = -3A<Ma 2 - cb 2 )\a 2 - <?\. 



(23) 
(x), obey- 



(24) 



The solution with finite energy obeys the BPS first order 
equations 



/,2\3 



(25) 



and the model topological solutions are then given by 



<P(x) 



s Va 



VI + Aa 4 a 



with a corresponding energy density 

a 6 A 



p(x) 



(26) 



(27) 



(1 + Aa 4 x 2 ) 3 ' 
The total solution energy may then be calculated as 

E = la 4 VXn. (28) 
8 

Defining the solution thickness by S = 2{a 2 \f~\)~ 1 we 
obtain 



0(00) 2 



0.701, 



(29) 



validating this thickness definition. 

In Fig. [5] the profile of the kink solutions ([TBI and 
(f2"6"]) are plotted for comparison. It is possible to observe 
the slower asymptotic behavior of the new kink solution 
when compared with the <fi 4 kink solution. In fact, for 
large positive values of x, solution (|2"6")l may be expanded, 
originating the polynomial behavior 



4>{x) 



8\xJ '128 



(30) 



while, as it is known, the cf> 4 model exhibits an asymptotic 
exponential behavior 



(j) s (x) 



1 - 4e 



-x/S s 



+ 2e 



-8x/S s 



+ 



(31) 



Note that the slower asymptotic behavior of the new 
kink solution is an intrinsic characteristic of the model, 
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FIG. 2: The <f> 4 kink solution profile and respective energy 
density, solid line, and the new-kink (|26[) solution and energy 
density profile, dashed line, for a = b = A = A = 1. 



regardless of its solution thickness. The condition 

( 2) 

Ki \4>{ x —> i 00 )) = is associated with the asymp- 
totic polynomial behavior, which is always slower than 

the exponential one obtained for >(<f)(x -> ±oo)) ^ 0. 

The new-kink energy density, shown in Fig. [51 also 
exhibits a polynomial behavior 



p(x) 



a 6 A 
~6T 



while the energy density of the 
usual exponential behavior, 



, \x\ » 1, (32) 



solution presents the 



p s (x) = 166 4 Ae 



-8x/S s 



l-4e 



-4x/S s 



+ 



X\ > 1. 



(33) 

For completeness, in Fig. |31 the model solutions asso- 
ciated with non-vanishing pressure values are presented 
in the ((j),d(j)/dx) plane. 

It is possible to calculate the force between a kink 
and an anti-kink solution, spaced by a distance L much 
larger than their thickness, using the energy momentum 
density T 01 component The resulting force is given 
by |F| cx (L/5)~ 6 for the new model solution, with E 
given by Eq. (|2"51) . clearly more intensive than the usual 



Yukawa-like e~ 2L / Sa , obtained for (j) 4 model. The com- 
parison of the two type of kink interactions was confirmed 
by numeric simulating the full time dependent field equa- 
tion, given by Eq. ([5]) for the (fT4"]) and (|2H|) potentials, 
with X = X = a = b = l. The simulation starts with a 
well separated kink/ anti-kink pair at rest, for both the 
new-model and 4> A model solutions. The corresponding 
initial configuration is 0o = <f>(x + xq) — <f>(x — xq) — 1. 
The results are presented in the upper panel of Fig |H 
where snap-shots of the two model solutions are shown 
for the same time-steps. It is possible to see that the 
time-scale associated with the two forces is very differ- 
ent. At the simulation time-scale the new model pair 




FIG. 3: Solution in the (<j>,d<f>/dx) plane for the cj> 4 model 
(|14p (left panel) and the new model (|20|l (right panel). 
The solutions for positive pressure values are represented in 
blue (dash-dotted lines), for negative pressure values in red 
(dashed lines), and for a vanishing pressure the solution is 
represented in black (solid lines). 



of solutions starts collapsing while the corresponding <fi 4 
pair does not move significantly. The different time-scale 
of the two forces is even more clear in the lower panel of 
Fig. [H where the positions of the two anti-kink solutions 
for the new model (left) and for the </> 4 model (right), are 
plotted with time for the two models. 



IV. STABILITY 

In order to investigate behavior of the solution un- 
der small perturbations, the following general fluctua- 
tions of the scalar field are introduced: 4>(x, t) = (j)(x) + 
Si Vi( x ) cos(wi t), where r]i(x) represents a set of pertur- 
bations of the static solution. Introducing these fluctua- 
tions, expanded to the first order in rji, in the equation of 
motion, we obtain an equation of motion formally anal- 
ogous to the stationary Schrodingcr equation, dubbed 
Schrodinger-like equation: 



dx 2 



+ U(x))r ]l = Ufa Hrn = u 2 ^ (34) 



with 



U{x) = V^\ = ^ {x) 



(35) 



An eigenfunction corresponding to u> 2 > would be a 
localized shape change of the kink, characterizing the ex- 
citation of an internal mode [4lj |. The formal analogy 
with quantum mechanics will show that the kink acts as 
a potential well for the linear waves 41J, and allows us 
to further in this paper apply the probabilistic interpre- 
tation for the study of frequency modes and resonances. 
This interpretation is standard in problems of field local- 
ization in branes (see, for instance, Refs. [42|, EH). 
For models that support a superpotential function, 

,2„ 



equation (|3"4"]l can be written as S^Sij 
St = — ! \Y 



d 
dx 



(2) 



and S 



d 
dx 



cj^rji, where 

r(2) 



(36) 
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FIG. 4: Upper panel: Numerical simulations of full time 
dependent field equation of the kink/anti-kink for the new 
model (solid line) and for the 4> 4 model (dashed line). The 
higher intensity of force for the new model is related with a 
stronger and consequently faster kink/anti-kink interaction. 
Lower panel: Corresponding positions of the anti-kink center 
for the new model (left) and for the (f> 4 model (right) with 
time. 



The non-negative value of H = S^S implies that the 
system does not support any state with a non-negative 
eigenvalue, ensuring the linear stability of the model. In 
terms of the superpotential function, Eq. (|35[) may be 
written as 



U{x) = (W i 



(2)n2 



(37) 



Using Eqs. (|T5j) and (I57|) . the Schrodinger-like potential 
for the 6 A model can be obtained, 



U(x) = 2b 2 X (2 - 3scch 2 (&VXr) 



(38) 



This potential is a modified rcflectionless Poschl- Teller 
44], that allows, in addition to the zero mode, a fre- 
quency bound mode with eigenvalue 2/3 2 A. All the others 
states of the model, with wf > 4b 2 A, are unbound. 

Equivalently, it is possible to obtain the Schrodinger- 
like potential for the new model, using the Eqs. (j2"Tj) and 



(EZJ), 



U(x) = I2\a 



4 (a 2 VAx) 2 -l/4 
[l + (a 2 y/Xx) 



212 



(39) 



This potential, shown in Fig. [5j is called the Volcano 
potential. Note it goes to zero at infinity and exhibits 
two maxima points with the value U max = 12Aa 4 /5, for 
4> = ±VlEa/5. 





'U(x) 
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FIG. 5: Profile of the Volcano potential U(x) given by Eq. 
I|5gjl for a = A = 1. 



Using the superpotential function, Eq. (|2"Tj) , we obtain 
— 3\/A0 \f o? — (j> 2 which can be used to derive 



(2) _ 



Eq. ([3l)]) for the new model, 



3Aa x 



dx l + (a 2 VAx) 2 



(40) 



ensuring the non-existence of unstable bound modes with 
negative eigenvalues. Since the potential vanishes for 
large values of x, the only bound mode solution is then 
the zero mode: 



Vo(x) = 



5 7A 



(41) 



All the vibrational modes with k 2 > are asymptotically 
plane waves. The absence of positive bound modes for 
the new kink solution could in principle affect the rate 
of loss of energy by radiation in a kink-antikink collision 
process in comparison to the usual (f) 4 kink 25, [26| . How- 
ever, note that for eigenvalues with w 2 < U mlix a possible 
resonance structure could localize the energy density in- 
side the defect in a sufficiently long time to act as an 
effective bound state for the interacting properties of the 
kink. The influence of this effect for the formation of two- 
bounce windows (27| in a collision processes is currently 
being under consideration. 

The volcano-like potentials were also investigated in 
the context of the braneworld scenarios where gravity is 



6 



coupled to one real scalar field. For some models, the 
behavior of the massive gravitational modes inside the 
defect are described in 1451-44811 . 



V. GENERALIZATION 

A generalization of the previous model, for solutions 
with a progressively slower power-law asymptotic behav- 
ior, is given by 



V 



A 



,2|n+l 



(42) 



where n is an integer parameter. A similar extension for 
the model of Ref. [4y| was considered in Ref. [49| in 
the context of the deformation theory. We see that the 
generalized potential preserves the maximum and min- 
ima values for <p max = and 4> m in = ±a. The first non- 
vanishing derivative orders of the potential at the minima 
are given by |U ( " +1) (±a)| = 2™ (n+ l)!Aa" +1 . Note that 
the cases n = 1 and n — 2 correspond to models (|T4"]) and 
(|2U)l respectively. The potential is plotted in Fig 6(a) for 
different values of n. It is possible to observe the enlarge- 
ment of the minima flatness with the increasing value of 
n. 

The corresponding generalized superpotential function 
W(4>), for \cj)\ < a is given by 



W{4>) = \a n ^<)>F 1 -,--( n + l); 



1 1 



2' a 2 



for \cj)\ < a. 
(43) 

Using this potential we may calculate the total energy of 
the generalized solution as, 



E,, 



: r (f 



f) 



r(2 + |)- 

Considering the generalized equation of motion, 



dx 2 



= -(n + I)A0(a 2 - 4> 2 )\a 2 - <\> 2 \ 



(44) 



(45) 



we determine the numeric solutions of the generalized 
class of models, and their corresponding energy densities, 
plotted in Fig |6(b)| and 6(c) for different values of n. Is is 
clear that the solutions and energy densities asymptotic 
behavior slow down with an increasing n. 

Using a perturbative method we also calculated the 
analytic asymptotic behavior of the generalized model 
solutions and corresponding energy densities, for different 
values of n > 1, 



4> n (x) = a 



1 - 



2x, 



(46a) 



2 <"+ 1 > / S n \ 



Pn (x) = a' n+ \n - !)- — [ f\ +... (46b) 
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FIG. 6: Higly interactive kink generalized class of models, 
Eq. H42|) . The dashed-dotted and dashed lines refer to n = 1 
and n — 2, respectively. The solid lines refer to n = 3 to 
n = 9. We have chosen A = a = 1. 



where we define S n = 2(^/~\a n )~ 1 as the n th solution 
thickness. Note that the asymptotic behavior of these 
solutions is included in power-law behavior presented in 
Eq. © at section HI 

For a better illustration of the progressively slower 
asymptotic behavior of the solutions we have also plotted 
the ratio between the value of (f> for x = S n /2 and its value 
at the minima, 4> n {5 n /2) / 4> n (oo), in Fig. [7] (asterisks). 
A loss of accuracy in the thickness definition, <5„, could 
be associated with an increasing of the solutions spread, 
however this behavior is accompanied by an increasing 
localization of their energy density, as shown in Fig |6(b)| 
This result is also confirmed in Fig J7] where the solutions 
energy densities (circ les) are plotted for x = S n /2, and 
moreover in Fig 6(c) (boxes) where the ratio 



=S n /2 



=-5„/2 



p n (x)dx 



(47) 



p n (x)dx 



x — ~oo 



is shown, corroborating the increasing of the energy den- 
sity localization with n. 

Finally in Fig. |6(d)| the Schrodinger-like potentials 
U n (x) are plotted for different values of n. The Vol- 
cano potential well tapers with the increase of n and 

n 



the maxima points in 
2"- 1 n(n - 2)™- 2 (2n - I) 2 -"Aa™. 
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FIG. 7: The ratio (j> n (5 n /2)/<j> n (oo) (asterisks), the energy 
density of the solutions p n (5n/2) (circles) for x — S n /2 and 
the energy ratio given by Eq. (I47|) (boxes), for different values 
of n 
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FIG. 8: Single resonance peaks characterizing the relative 
probability P re i{Q) for finding frequency modes near to the 
center of the kink, as a function of the frequency to of the 
mode. The peaks are for n — 2 to n = 9. The lower is n, the 
thinner is the corresponding resonance peak. 



with n may be associated with occurrence of resonances 
with higher frequencies when n is increased. In order to 
confirm this hypothesis, we considered even parity states 
and applied the well-known Numerov method [50| (for 
application of the method to brane problems see ref. [471 ] 
and references therein) for finding the frequency modes 
r] m (x). We want to compare the relative probability 
P re i(0) for finding a frequency mode around the kink 
center located at x = 0. This justifies choosing a fixed 
control distance x c and defining, for a known numerical 
Schrodinger-like potential for 



Prez(0) 



r] m {x)dx 



(48) 



r) m {x)dx 



where x max allows for a sufficiently large number of oscil- 
lations in the range x rna x < x < x max for the wavefunc- 
tions 7] m (x). Fig. [8] shows the resonance peaks P re i(0) 
found for n — 2 to n — 9, as a function of the frequency 
uj of the mode. The first thin peak is for n — 2. The 
sequence of peaks in the figure follow the increasing se- 
quence of n. Note that each peak characterizes a res- 
onance centred at a particular value of m. The whole 
sequence shows that the increasing of n turns the reso- 
nance peak broader. The first thinner and higher peak 
for n = 3 is the first modification to the 4 model that 
has a bound frequency state. The modification for n = 3 
trades the bound state for n — 2 by a quasi-normal state 
with lifetime inversely proportional to the peak thick- 
ness. The growth of n leads to higher frequency resonant 
states, agreeing with our discussion after Eq. (|47p . 



VI. ENDING COMMENTS 

In this work we have introduced a new real scalar field 
model that admits kink solutions exhibiting a power-law 
asymptotic behavior, instead of the usual exponential 
one. This is a consequence of having a potential with 
a null second order derivative at the minima. We have 
shown that this particular asymptotic behavior of the 
solution leads to a much stronger kink/ anti-kink inter- 
action, when the same thickness and energy-scale of the 
defects is considered, than in the standard models. 

We have further generalized the model to a class of 
models admitting n vanishing derivative orders of the 
potential at the minima. It was found that, despite the 
increasing of the solutions spread with n, their corre- 
sponding energy density becomes more localized when n 
is incremented. 

Stability analysis led to a Schrodinger-like equation 
where bound modes are absent. This may be related 
with interesting properties of defect collisions. The two- 
bounce effect [25|, |27| appearing in the collision of (j) 4 
model defects requires the transference of a portion of 
kinetic collision energy to the internal positive bound 
mode of the solutions. Moreover we have shown that the 
Schrodinger-like potential of our model is the Volcano 
potential. The inexistence of frequency bound modes in 
this potential suggests that no two-bounce effect is ob- 
served in this case, however a resonance effect may occur 
inside the defects during the collision, leading to some 
peculiar behavior. This is an interesting future work. 

Moreover the embedding of these kink solutions in four 
spacetime dimensions generates a domain wall network 
with a particular dynamics, which could have interest- 
ing applications in Cosmology or in Condensed Matter 
in the non-relativistic limit. Note that different energy 
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factors have to be considered in dynamics of domain wall 
networks, not only the kinetic energy or the radiation 
energy losses but also the wall curvature needs to be put 
into play. Possible implications of this class of models in 
the dynamics of domain wall networks is also an inter- 
esting future study. 

Another possible relevant application for this class 
of models is in the braneworld scenario. Indeed, as 
one knows, using the first-order formalism the exten- 
sion for (4, 1) dimensions of this kink-like solutions is 
straightforward (l6j |. Gravity and particle localization 
in branes have been studied extensively in the literature, 
with the main focus of looking for resonance phenom- 
ena [H |H, m liHHsl. For the case presented here, the 



special characteristics of perturbations around the kink 
solution suggests that similar studies is also a very inter- 
esting issue to be further considered. 
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